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i ^ ! The particle production in the self-interacting iV-component complex scalar field theory is studied 
£,t large N. A non-Markovian source term that includes all higher order back-reaction and collision 

t-h effects is derived. The kinetic amplitudes accounting for the change in the particle number density 
caused by collisions are obtained. It is shown that the production of particles is symmetric in the 
>■ momentum space. The problem of renormalization is briefly discussed. 

m ; 

O : PACS: ll.lO.z; 05.20.Dd; ll.10.Ef; 98.80.Cq 
O 

o 



Introduction 



Qjn recent years an essential progress has been achieved in establishing a consistent kinetic description 
J-h Of particle production in quantum field theory from first principles. Investigations performed for 
^ QED with a constant [1] or time-dependent external electric field [2, 3] have shown that a consistent 
• '-j field-theoretical approach leads to a kinetic equation with a modified source term, providing a non- 

^> Markovian evolution of the distribution function. 

H , 

Non-Markovian effects in collision terms have been studied in connection with relativistic heavy 
ion collisions [4], collective effects in nuclear matter [5], nuclear fragmentation [6], and the damping 
rates of giant dipole resonances [7]. The kinetic description of transport phenomena in QCD [8, 9] 
and scalar field theories [10] has been also developed. A method to obtain quantum kinetic equations 
directly from the underlying quantum field theory implementing the Schwinger-Keldysh formalism 
and a dynamical renormalization group resummation has been presented in [11] . This method has 
been applied to different scalar and gauge theories, including the self-interacting <y9 4 -theory, yielding, 
however, no memory effects. 

In the present paper, we aim to derive a non-Markovian quantum kinetic equation for a self- 
interacting scalar field theory. We consider the model of a massive iV-component complex scalar 
field in the large- N approximation. Our purpose is to obtain a non-Markovian source term which 
incorporates all higher order, i.e. back-reaction and collisions effects. 

In contrast with [11] , in our paper we follow a kinetic approach introduced recently to study the 
decay of CP-odd metastable states and based originally on the evolution operator technique [12] . 
Herein we generalize this approach to include higher order effects. 

1 e-mail: fuad_saradzhev@hotmail.com 



1 



Our plan is as follows. In Section 2 we introduce the model and perform the quantization of 
the scalar field fluctuations within the mean-field approximation. We construct the fluctuations 
Hamiltonian and identify its diagonalized part. In Section 3 we study at large N the effects of 
the remaining nondiagonalized interactions in the particle production. We derive a quantum ki- 
netic equation governing the production process and discuss its renormalization. We conclude with 
summary in Section 4. 



2 Quantization 

The Lagrangian density of our model is 

N 



£ = E [V&Wvi) - - v(<p), (l) 

where m is a bare mass. The model is defined in a finite volume L 3 , —L/2 < X{ < L/2, i = 1,2,3, 
and the scalar field fulfils the periodic boundary conditions. In the continuum limit J2k > where 

the summation is over discrete momenta k = ^-n, (ft) = (ni,n 2 ,n 3 ), goes to J -^s- 
The potential is taken as 

^)=2(E^Vi)» A>0 > ( 2 ) 

3 

and includes self-interaction for each component as well as interaction of different components. From 
(1) we obtain the Klein-Gordon type equation of motion 

{n + m 2 )<p j = J j (3) 

with the current 

dV N 

J j = -^t = -*(X>iViW ( 4 ) 

9<Pj 1=1 



2.1 Mean- field approximation 

Following the mean- field approximation, we usually decompose ipj(x,t) into its vacuum mean value 
4>j(t) = (ipj(x,t)) and fluctuations Xj 

<Pj(x,t) = </>j(t)+Xj(x,t) (5) 

with (xj(x, t)) = 0. The mean field is treated as a classical background field defined with respect to 
the in- vacuum |0) as 

Ut) = js /^<0|^(£,*)|0>, (6) 

so that in the limit t — > — oo (j>j(t) — > 0, while the fluctuations are quantized and take place at all 
times. 

The background 0-fields are in general complex and contribute to the charge of the system. By a 
global U(l) gauge transformation the phases of the 0-fields can be put equal to zero at a fixed time. 
If the vanishing of the phases is compatible with equations of motion, it is possible to make them 
zero at all times as well. Otherwise, the background phases become dynamical variables coupled 
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to the fluctuations in all orders. This results in new, non-diagonal terms in the Hamiltonian and 
therefore in new interactions. 

To simplify the problem, it is convenient to choose the background neutral. Let us introduce the 
real fields 

<^ = ^- + 0*), (7) 
1 

2i 

and rewrite the decomposition (5) as 



2j = vW-tj)' 



(fj (x, t) = y (t) + Xj (x, t), (9) 

where Xj = Xj+ifoj- The vacuum mean value of the "new" fluctuations is non-vanishing, (% ■) = i&j- 
Using Eq.(9) provides the following decomposition for the current 

N 

jj ( <p) = J j (0i ) - nlxj - i*j w (xi + x\ ) + Jj (0; x) , ( 1 °) 

i=i 

where 

= ( n ) 

^- = >/A0ij, (12) 
while Jj((p;x) includes terms of second and higher orders in fluctuations, 

1^-x) = -VX^xki - V^XjEMXi + x]) + MX)- (13) 
i=i i=i 

Substituting Eq.(9) also into (3) and taking the mean value (...) yields the vacuum mean field 
equations 

fa, + (m 2 + fiDfa = Re(J i (0; x)>, (14) 

4>2j + ijn 2 + /4)(f>2j = lm(Jj((f);x)}- (15) 

As seen from (15), only if higher orders are neglected can the 0-fields be taken neutral. Assuming 
that (f> 2 j = leads to the condition 

Im(7 J (0;x))(02 = O) = O, (16) 

which is not valid in general if the fluctuations are charged. 

Eq.(3) in concert with Eqs.(14)-(15) provides also the equations of motion for the fluctuations. 
Introducing instead of Xjix] the hermitian fields 

Xij = ^(Xj + x]\ ( 17 ) 
1 



with (xij) = 0, (x 2 j) = <hj, y ields 



X2 3 = —j^{Xj-X\), (18) 
{o + m 2 aj )x aj = F aj , (19) 



(a = 1,2), where 



are effective masses, and 



mfj = m 2 + /2 2 + 2/4 (20) 

"<lj 



m 2 = m 2 + i2 2 (21) 



N 

F lj = 7 y - (7i,> - 2 N Yl ViXw (22) 

i=i 
¥j 

F 2j = J 2 j, (23) 

the hermitian currents J a j being defined instead of Jj, j] in the way analogous to (17)-(18). 

Since the effective masses squared are always positive, the system evolves only in the non-tachyonic 
regime ( for the kinetic approach to the tachyonic regime, see [13] ). For N — 1, the effective masses 
squared take the form m\ = m 2 + 3A0^, m 2 , = m 2 + \(f> 2 in agreement with [14]. 

Rewriting (19) for the Fourier components x a j(k, t), a = 1, 2, we obtain a Mathieu type equation 

t aj (k,t) + u% (k, t)x aj (k,t) = F aj (k, t) , (24) 
where F aj (k,t) is the Fourier transform of F aj (x,t), and 

u; 2 aj = k 2 + m 2 aj (25) 
are the time-dependent frequencies squared of the fluctuations. In the in-limit, u a j become 



lim w Qj (fe,t) = uj°(k) = \j(k) 2 + m 2 . (26) 

Eqs. (14), (15) and (24) are self-consistently coupled and include all higher-order effects. The vacuum 
mean fields modify the equations for the fluctuations via the time dependent frequencies, while the 
fluctuations react back on the 0-fields via the source terms Re(Jj) and Im(Jj) in Eqs.(14),(15) and 
on the fluctuations themselves via the term F a j(k,t) in Eq.(24). 



2.2 Hamiltonian 



The fluctuations act on each other in two ways: directly, when higher orders are included, and via 
the 0-fields. The potential (2) is decomposed as 

N 2 N N 

v(<p) = vifa) + yftfi £ than + v + xh) + E vmxijxu + v(fr x), (27) 

j=i i=i (j,0=i 

where V(4>; x) contains third and fourth orders in fluctuations, 

N 

v{bx) = - E^(4(x) + J}{X)) + V(x). (28) 
j'=i 

With the decomposition (27), we deduce from (1) the Hamiltonian density governing the dynamics 
of the fluctuations 

Y 2 N N 
K x = o E E «• + (Vx aj ) 2 + m^x^) + E WiXijXu + «x> (29) 

a=lj=l (j,0=i 
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where the momenta canonically conjugate to x a j-> anc ^ 

N 

n x = ¥(<!>,■, x) + J2(J^)Xir (30) 
J'=l 

The expression (29) is not diagonal even in the second order in fluctuations. The non-diagonal terms 
include interactions of fields of different values of j. 

In terms of the Fourier components % aj - W a j(k,t), the fluctuations Hamiltonian takes the 
form 

N N _ 

H x = / d 3 xH x = E ^ + E H Ji + L 3/2 H x (k = 0, t), (31) 
i=i 0,0=1 



where H x (k,t) is the Fourier transform of the potential H x (x,t), and 

2 

2 



H i ■> = \ E E (<*(M)^-(M) + <(MxUMx Qj (M)), (32) 



a=1 jfc 



= W E OxW^ 0- (33) 

The Hamiltonian equations of motion are 



tajik, t) = t4-(M), (34) 
# aj (M) = -^.(fc,*)x^-A0 + ^i-A*)> (35) 

in agreement with (24). 

By the transition to the new fields (17)-(18), we have achieved a partial diagonalization of the 
Hamiltonian in the second order in fluctuations. Its complete diagonalization, including all higher 
order, i.e. interaction terms, is a very complicated problem. However, there are different approxima- 
tion schemes which allow one to diagonalize a "physically" important part of the interaction, while 
the remaining part is treated as a small perturbation. 

Let us assume that, by using one of these schemes, we succeeded in diagonalising some of the 
interaction terms in our model. We assume next that a new diagonalized part of the Hamiltonian 
looks like the one in Eq.(32), the frequencies of fluctuations being only modified. Namely, u^(k,t) 
is replaced by 

^.ft0 = <fi0+Mj(^,0, (36) 

where M^(k, t) represents the higher order effects contribution. The explicit expression for M^(k, t) 
depends on the approximation scheme under consideration. In the next section, we will calculate 
M^(k,t) in the large- N approximation. 

For the Hermitian fluctuations, we introduce then the annihilation and creation operators d a j(k, t) 

j. — * 

and d a j(k,t) in a standard way by 

X aj (k, t) = T aj (k, t)d aj (k, t) + F* aj (k, t)dl 3 (-k, t) (37) 

and 

W aj (k,t) = -iujaj(k,t) T aj (k,t)d aj (-k,t) - r* aj {k,t)d) aj {k,t) , (38) 
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where 

r aj (k,t) = = = :exp{-iQ aj (oj aj ,t)} (39) 
\j2uj aj (k,t) 

and Q a j(uj a j,t) are arbitrary phases. In the case when higher order effects are omitted , i.e. M\ - = 0, 

the phases Q a j(u a j,t) take the form u k t in the in-limit, and the operators Xaj(k,t), ir a j(k,t) can 
be connected with the corresponding in-field ones by making use of the evolution operator [12]. 

The ansatz (37) ( and (38) ) has the same structure as the free field theory one except the 
frequencies u a j(k,t) replace u^(k), and the time- independent annihilation and creation operators 
are replaced here by the time-dependent ones d a j(k,t), d^{k,t). The Hamiltonians Hj become 

1 2 

H i = 9 EE^'(M)(i(M)^(M) + <MM)4*(M))- (4°) 
3 Kinetic equations 

There are 2N types of neutral particles in our model which are related to 2N types of the fluctuations 

— * 

X a j- The numbers of these particles of momentum k at time t are given by the occupation number 
densities 

a/;,(m) = (oi^mk^OIo)- (4i) 

In the limit t — > — oo, these densities vanish because the initial state is assumed to be "empty", i.e. 
without particles. 

In general, N a j{k,t) and J\f a j(—k,t) are not equal. Therefore it is convenient to introduce 

Af aj ,±(k, t) ee ^(Af aj (k, t) ± Af aj (-k, t)), (42) 

— * — * — * — * 

where J\f a j >+ (k,t) is the particle number averaged over the directions k and (—k), while J\ a j_(k,t) 
measures the degree of asymmetry. 

We consider first the time evolution of J\f a j-(k, t). Using the relations 

2r ai (/c,t) v v aj (k,t) 

and the Hamiltonian equations of motion (34) — (35), we obtain 

A^,-(M) = -Im(0|x^(fc,t)F a ,(fc,t)|0). (44) 
If we omit higher-order effects and take J a j = 0, then 

N 

AT aj ,-(k,t) = 2 f x j 6 al Y,ViMQ\Xi j (k,t)x ll (k,t)\0), (45) 

i=i 

i.e. the time evolution of the asymmetry in the production of particles by the field of a fixed value of 
j is determined by its interaction with other fields. For N — 1, the density H a \,-{k,t) is conserved, 
Jv a i-(k,t) = 0, so if we take J\f a i,-(k,to) = 0, where t is a moment of time at which the particle 
production starts, then the asymmetry does not appear at any time t later. 
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Taking next the time derivative of J\f a j !+ (k,t), yields the evolution equation 



Afcy,+ (M) 



CO 



Re 



e- 2ie ^C aj (k,t) + A%(k,t), 



ma 

aj 



where we have defined the time-dependent pair correlation functions 

C aj (k,t) = (0\d aj (-k,t)d aj (k,t)\0) 

and 



—Re(0\lT aj (k,t)F aj (k,t)\0) 



uj, 



with 



(46) 

(47) 
(48) 
(49) 

, (50) 
(51) 

As seen from Eqs.(46) and (50), the higher order effects contribute to the evolution equations via 



F aj (k,t) = F aj (k,t) + M 2 aj (k,t)x aj (k,t). 



The pair correlation functions C a j(k,t) can be shown to obey the equations 
C aj (k, t)-2i(O aj -u aj )C aj {k, t) = [^L(l+2N aj , + $, *)) + (- Ag(fc, t) +tA$(k, t 



where 



Ag } (M) = Re(0|xL,(M)i^(M)|0>. 



the frequencies oj a j and via the functions A^(k,t) and A^(k, t). These functions include higher 
order correlation functions, namely, the vacuum expectation values of the products of three and four 
field operators. 



(2), 



3.1 Large- N approximation 

To calculate A^(k,t), A^(k,t), we use the large- N approximation [15-18]. At large N, the theory , 
which consists of an infinite number of coupled correlation function equations, reduces to the solution 
of coupled equations for the one- and two-point functions. For simplicity, we consider the case of 
spatially homogeneous vacuum configurations and assume that (0|x 2j -(^, jO) is a function only of 
time, then 

<O|x 2j O?,t)|O> = 2r (52) 

The basic points of the large iV-approximation are as follows, (i) The equal time two-point correlation 
functions are proportional to 5 a ^ (a, (3 = 1,2) and 5ji = 1,N) unless the vacuum expectation 
values of x a j are non-zero. Assuming translational invariance of the vacuum state |0), we define the 
correlation functions 

(0\Xij&t)Xu(M\V) = ^G l3 (x-y;t), (53) 
(0\x 2j (x,t)x2i(y,t)\0) = <j) 2j <p2i + SflGzjix-fit), (54) 

and 

(0|7f y (f,t)Xii(j7,0|0) = SjiD^x-fit), (55) 
(0\w 2j (x,t)x 2 i(y,t)\0} = <j) 2 jhi + 8jiD 2j (x-y;t), (56) 

while 

(0\Xi 3 (x,t)x 2l (y,t)\0) = (0|7fy(f,t)x 2 i(j/,*)|0) = (0|7f 2j (f,0Xi/(y,i)|0) = 0. (57) 
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(ii) The vacuum values of all fields are the same, i.e. (fiij = 0i and <p2j = 4>2 for all values of j. 
As a consequence, /z 2 , = A0 2 where A = AiV, and the effective masses coincide, m\j = m\j = m 2 = 
m 2 + A0 2 . The correlation functions G a j, D a j as well do not depend on the particular value of j. 
However, in what follows we keep the index j on all correlation functions, vacuum mean values of 
currents and etc, except the 0-fields. 

(iii) The higher order correlation functions are factorized into products of two-point correlation 
functions and vacuum fields. For example, 

(0\X2j&t)Xii(x,t)Xu(x,t)\0) = <hG u (P;t), (58) 

(0\x 2 j(^t)x 2l (x,t)x 2l (x,t)\0) = &[G 2 ,(0;f) + 2fyG 2j -(0;f) + 3$]. (59) 

With these points it is straightforward to calculate A^j(k,t) and A^j(k,t). We start, however, 
with the vacuum mean value of the hermitian currents J a j(x,t). Using (13) yields for these currents 
the following expressions 



N I N 

J y = -V2\xij E ViXu - 2 + A ^') E {xii + Xli) , (60) 



A A 



2' 

i=i z i=i 



J 2j = -V2X X 2j E ViXii - 1^X23 E (xu + xli) • (61) 
1= 

Using next the factorization (iii), we obtain 

(7 y > = -^n^i, (62) 



2 L" J 

where 

1 2 



(J 2 ^ = -^ 2 + 30 2 U 2 , (63) 



nj = G y (0;f) +G 2J (0;0 = — £ E<°lwM)x*(M)|0>. (64) 



As seen from Eq.(63), (J 2 j) = for 2 = 0, i.e. the condition (16) is fulfilled. We come therefore 
to the important conclusion: in the large- N approximation the 0-fields can be taken neutral at all 
times. Henceforth we put 2 = 0. 

Substituting Eq.(62) into Eq.(14), we rewrite the vacuum mean field equation as 

0i J + K 2 + ^?)^ = O, (65) 
so that the higher order effects contribution to the effective masses is 



Mf = (66) 



We turn now to the calculation of A^j(k,t) and A^j(k,t). In the large- N approximation, due 
to the point (i) the interaction of components with different values of j does not contribute to the 
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vacuum expectation values, and A a Uk,t) and A a j(k,t) take the form 



aS(M) 



U). 



-Re(0\w aj (k,t)J aj (k,t)\0) 



A 



2u a f 1 



Q 2 -ReD aj (k,t), 



A 



Ag(M) = Re(0|£,(^t)^(^t)|0) + ^ J 2 -G a ,(A;,t) ) 



(67) 
(68) 



where D a j(k, t) and G a j(k, t) are the Fourier transforms of D a j(x—y; t) and G a j(x—y; t), respectively. 
Evaluating next the self-interaction contribution by making use of (iii), we obtain 



Re{0\W aj (k,t)J aj (k,t)\0) = --Q 2 -ReD aj (k,t), 



A. 



Re(0\xUk,t)J aj (k,t)\0) = -^ 2 -G aj (k,t), 



2' 3 

2 n > 



leading finally to 



A<g(M) = Ajg(M) = 0. 



(2), 



(69) 
(70) 

(71) 



In the same way we calculate the vacuum expectation value in the r.h.s. of (44) and obtain that 

A^-_(M)=0, (72) 

i.e. in the large-iV approximation the particle production is symmetric in the momentum space for 
all times. 



With A 



(1) _ A ( 2 ) 



0, Eq.(50) is solved by 



C aj (k,t) = e 2ie - f dt'^(k,t')(l + 2Af aj , + (k,t')) 



2i(e aj (k,t')-e aj (k,t)) 



where 



G aj (k,t) = f dt'u aj (k,t'). 
Substituting (73) into (46), we obtain the following closed equation for J\f ajt+ (k,t): 



(73) 



(74) 



Af aj , + (k, t) = ^W aj (k, t) j 1 dt'W aj (k, t') • (l + 2Af aj , + (k, *')) cos [x aj (k- 1', t)] , (75) 



where 



W = ^ 



CO, 



are kinetic amplitudes, while 



x a j(k',t : t) — 2 ® a j(k,t) Q a j(k,t) 



(76) 



(77) 



is the phase difference. 

Eq.(75) is a complete quantum kinetic equation in the large- N approximation with all higher 
order effects included. It determines the time evolution of the number of particles of a fixed j and a 
fixed momentum k. 
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The background 0-fields are given by solution of the non-linear equation (65), while Vt 2 - must be 
determined self-consistently from 

1 2 



75 EE 



with 



R a j (k, t) 



L 3 

1 

' 2 



2u, 



{l + 2Af aj , + (k,t))+2R aj (k,t) 



to 



dt'W aj (k, t') (l + 2J\f aj:+ (k, *')) cos [x aj (k; t\ t) 



(78) 



(79) 



We have therefore a system of three coupled equations. 

The kinetic equation (75) has the following important features: (i) The source term, i.e. its r.h.s. 
is non-Markovian. Therefore calculating the particle density at any given instant requires a complete 
knowledge of the history of the production process. In addition, the integrand is a non-local function 

(k; t', t) ; (ii) The background 



x 



a j \ 



of time, which apparent in the coherent phase oscillation term cos 
0-fields do not contribute to the kinetic equation directly, but only via the frequencies of the quantum 
fluctuations, as evident in Eqs.(76) and (77); (iii) The production rate is affected by the produced 
particles' statistics, as seen in the statistical factor (1 + 2Af ajt+ (k,t)). In the case of fermions, this 
factor would be (1 — 2J\f ajt+ (k,t)). 

The kinetic amplitudes W aj (k,t) are decomposed as 



w aj = w aj + w r 



n 

aj i 



where 



and 



W aj = 



o.7 



<*3 



w n = 

«J ~ 



A 

2< 



\pl 
2< 



:)" 



(80) 



(81) 



(82) 



The amplitudes W a j(k,t) account for the increase in the particle number caused by the transition 
of energy from the background 0-fields to fluctuations. As soon as the energy transition stops, the 
amplitudes W a j(k,t) vanish. If, for example, the time evolution of the 0-fields starts at t — t and 
ends at t — t±, then W a j(k, t) = for all t > t 1 . 

In contrast with W a j(k,t), the new amplitudes W^(k,t) account for the change in the particle 
number density caused by the interaction of the particles produced, i.e. by their collisions, and do 
not vanish when the particle production by the background fields stops. The fluctuations interact 
with each other in two ways, directly and via the 0-fields. For t > t± , the time-dependence of the 
particle number density is only governed by collisions. 

Collisions also modify the expression for the dynamical phase. If collision effects are omitted, Q a j 
reduces to the adiabatic phase 



eg(M) = / dt'cu aj (k,t') 



to 



that takes the form uj^\t — to) in the in-limit. With collisions taken into account, the dynamical 
phase is no longer the adiabatic one and the connection with the in-picture no longer exists. 

The non-Markovian kinetic equation in the collisionless limit can be simply obtained from (75) 
by replacing W a j with W a j and Q a j with 0^. The kinetic equations in both cases, i.e. with and 
without collisions, have therefore the same non-Markovian structure. This is a result of the large- iV 
approximation. However, collisions introduce non-Markovian effects in addition to those already 
present in the case without collisions. Even in the "collision" regime for t > ti, the production rate 
is affected by the background fields induced particle production process from t to t 1 . 
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3.2 Renormalization 

All the divergences in our approach come from the expression for uj a j(t): 



u 2 



A 



(84) 



the sum over k in the definition of Vt 2 [Eq.(78) having both quadratic and logarithmic divergences 
which have to be absorbed by mass and coupling constant renormalization [18] . The subtraction of 
Sl 2 At Q ) from VL 2 (t) removes the quadratic divergence. With J\f a j j+ (k,t ) = R a j(k,t ) = 0, we rewrite 

(85) 

k 2 + m 2 R , (86) 



where 



m 2 R being the renormalized mass squared, 



R 



m 2 R = m 2 + \(t)\{t ) + ^fi-(to), 



A. 



and 



S s (t) = \ <fi(t) - <fi(t ) 



+ 



A 



n?(f)-n?(*o 



(87) 



(88) 



S 3 (t ) = 0. 

The remaining logarithmic divergence is removed by the coupling constant renormalization. The 
bare coupling constant and the renormalized one are related by the geometric series 



A = 



A 



R 



i-x R sy 



where 



Multiplying both sides of (88) by (1 — Xr5\), we obtain the finite equation 

X R 



S j (t) = X R #(t)-#(f ) 



+ 



n 2 (t)-n 2 (t ) + 2S\-s j (t) 



If we expand fl 2 (t) and Sj(t) in a power series in Xr, 

oo 

^(*)=^(0)(*) + E^' n ?(n)(*), 
n=l 



-Sj'C*) = H ^ ' S j(n)(t), 
n=l 

and choose the initial conditions for fi|(£) as 

fl 2 (t ) = f^ (0) (t ), 

fi ?(fe)( t o) = for k ^ 1 ^ 



(89) 
(90) 

(91) 

(92) 
(93) 

(94) 
(95) 
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then (91) reduces to the following set of equations 



S,-(2)(0 = ^(i)(*) + *A.fi,- (1) (0, 

^■(3)(0 = ^(2)(*)+^-% 2 )(0>-> 




(96) 



i.e. in the large-iV approximation we need only to calculate in the second order to identify the 
counterterms [18]. 

The renormalized vacuum mean field equation takes the form 



Expanding the amplitudes W a j and the dynamical phases Q a j in power series in Xr, it is also 
possible to deduce from Eqs.(75) and (78)-(79) a renormalized kinetic equation in a fixed order 
of perturbations as well as to develop a perturbation theory technique for the particle production 
problem. This work is in progress and will be reported elsewhere. 



1. For the model of a self-interacting complex ^-component scalar field, we have derived in the large- 
N approximation a quantum kinetic equation with a non-Markovian source term. This equation 
determines the momentum distribution of the standard, non-tachyonic modes produced in quantum 
fluctuations of the scalar field around its vacuum mean value. 

Our equation is complete in the sense that it includes all orders of the back-reaction and collision 
effects. Its form is compact and appropriate for the numerical study. We have obtained an explicit 
expression for the kinetic amplitudes related to collisions. These amplitudes account for the change 
in the particle number density in collisions in both "energy transition" and "collision" regimes. 

Our method of derivation of kinetic equations is not restricted to the model discussed. It can be 
applied to any self-interacting field theory, providing the approximation scheme to treat non-diagonal 
interactions is specified. We have shown that in the large- N approximation the functions A^- , 
vanish, so that the source term with higher order effects included has the same non-Markovian 
structure as the one with higher order effects omitted. However, beyond the large- N approximation, 
the functions A^- , A^ are expected to modify the structure of the source term as well. 

2. For the particle production induced by an external electric field, the non-Markovian effects are 
known to be important when the field is strong [3] . For weak fields, non-Markovian effects disappear, 
and the Markovian limit of the kinetic equation, which is defined by the neglect of memory effects in 
the source term, can be successfully employed. If the external field is switched off, then the particle 
production stops. 

For the self-interacting fields, the situation is completely different. The self-interaction takes 
place at all times and at all values of fields. It introduces memory effects, which appear even at lower 
orders of perturbations and can not be therefore neglected in any sensible limit or regime. Once 
the self-interaction is taken into consideration, the non-Markovian character of the source term is 
unavoidable. Even if the transition of energy from the background fields to fluctuations stops, the 
particle production process continues, being influenced essentially by the background field evolution 
history. 




(97) 



4 Discussion 
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